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I.  INTRODUCTION 


In  the  past  several  years  the  US  Army  Missile  Command  (MICOM)  has  been 
engaqed  in  the  development  of  techniques  for  improving  the  self-alignment 
performance  of  an  inertial  measurement  unit  (IMU)  [l  to  12  J.  This  report 
documents  the  continued  findings  of  new  techniques  applicable  to  IMU  self¬ 
alignment. 

The  precision  of  and  the  required  time  for  IMU  self-alignment  depend 
not  only  on  the  quality  of  the  IMU  hardware  but  also  on  the  effectiveness 
of  the  attendant  software.  In  fact,  an  ingeniously  developed  self-alignment 
algorithm  enables  the  compensation  of  anomalous  effects  caused  by  hardware 
imperfections.  Thus,  the  stability  of  the  hardware  characteristics  has 
become  more  important  than  the  exactness  of  the  hardware. 

In  general,  an  IMU  self-alignment  procedure  contains  on-line  determi¬ 
nation  of  a  set  of  parameters  characterizing  the  state  of  the  IMU.  To  cope 
with  the  ever  existing  effect  of  random  noise  and  disturbances,  redundant 
measurements  and  a  regression  type  of  data  reduction  technique  are  used  for 
improving  parameter  determination  accuracy.  The  purpose  of  this  report  is 
to  describe  two  features  of  data  reduction  for  the  self-al iqnment  of  a  gim- 
baled  IMU.  The  procedure  involves  an  accelerometer  based  two-position  self- 
aliqnment  technique  which  has  been  reported  in  detail  elsewhere  [l,  5,  7]. 
The  two  features  to  be  discussed  are  a  continuous-time  algorithm  and  a 
"scraping"  technique.  Since  the  present  concern  is  the  data  reduction  soft¬ 
ware,  only  the  necessary  background  will  be  reviewed. 

II.  THE  ANALYTIC  MODEL 

The  parameters  to  be  determined  during  IMU  self-alignment  are  contained 
in  the  following  two  measurement  equations, 

VN(t)  =  Ajt  +  A2t2  +  uA4t3  (1) 

VE(t)  =  A3t  +  A4t2  -  uA2t3  (2) 

where  A.,  i=*1  to  4,  are  the  parameters,  t  is  the  time  variable,  P  is  a  con¬ 
stant,  and  Vf!(t)  and  VE(t)  are  observed  north  and  east  velocities.  Bv 
letting  t«kT,  the  discrete-time  countemarts  of  Eouations  (I)  an'1  (2)  are 
obtained  as 

vN(kT)  =  Bjk  +  B2k2  *  vB4k3  (3) 

vEfkT)  =  B3k  +  B4k2  -  vB2k3  (4) 
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where 


B1  =  AjT  B2  =  A2T2 

B3  =  A3T  B4  =  A4?2 

v  =  uT. 

The  discrete-time  model,  given  by  Equations  (3)  and  (4),  was  used  by 
the  present  authors  in  the  previous  development  of  a  data  reduction  algo¬ 
rithm. 

III.  CONTINUOUS-TIME  REGRESSION  ALGORITHM 

The  continuous-time  regression  algorithm  differs  from  the  discrete¬ 
time  one  in  that  the  former  uses  integrations  while  the  latter  uses  sum¬ 
mations.  It  is  apparent  that  the  latter  is  the  approximation  for  the  former. 
A  simple  demonstration  of  this  effect  is  given  below. 

Example:  Consider  a  polynomial  model 

V(t)  =  at  +  bt2  +  ct3.  (5) 

and  the  associated  integral  square  error 

l\  =/o  [V(t)  "  (at  +  bt2  +  ct3)]2dt.  (6) 

The  discrete-time  versions  of  Equations  (5)  and  (6)  are 

V ( k )  =  axk  +  br2k2  +  cr3k3  (7) 

and 

N 

U  =  l  t [v ( k )  -  (ark  +  bx2k2  +  ctV)]2  (8) 
k=l 

where  N  -  t/t.  Elements  of  the  matrices  of  the  associated  normal  equations 
are 

Cm  =  m  =  2  to  6  (9) 

for  the  continuous-time  model,  and 

c_  =  Tm+1  £  km,  m  =  2  to  6  (10) 

m  k=l 
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for  the  discrete-time  model.  Let  t=240  and  i=.192,  then  N~ 240/ . 192=1250. 
The  computed  values  of  Cffl  and  are  listed  in  Table  I.  Percent  errors  in 

Cm  caused  by  use  of  the  discrete-time  model  are  also  listed  in  the  table. 

Return  to  the  I MU  measurement  model  represented  by  Equations  (1)  and 
(2).  The  integral -square  error  resulting  from  a  continuous-time  least- 
square  regression  is  given  by 

UV  ’  Jk [yo  -  (*,t  *  A2t2  f  ,,A4t3}]2dt 

*iok[',E(«  -  (A3t  *  a/  -  „A2t3)]2dt.  (11) 

Letting  *T  =  41  ■  41  =  41  =  0 

i\]  lA,  4A3  4a4 

yields  the  following  set  of  normal  equations  which  specify  the  condition 
for  I(tk)  to  be  minimum,  where  all  integrations  are  from  t=0  to  t=tk. 


Hdt 

=  Aj/t2dt  + 

A2Jt3dt  + 

(12) 

J(,2vn 

-  ut3VE)dt  = 

=  Aj ^t3dt 

+A2(|t4dt  + 

q2/t6dt)  - 

A3u/t4dt. 

(13) 

JtVEdt 

=  A3Jt2dt  + 

Ajt3dt  - 

pA^t  dt. 

(14) 

|!t2vE 

*  ht3VN)dt  = 

=  A3Jt3dt 

+  A4  (Jt4dt 

+  u2/t6dt)  +  A1ujt4dt. 

(15) 

Uefi ne 

Wj  =  JtVt 

W2  -J(t\  -  ut3VE)dt 

W3  =JtVEdt  (16) 

W4  =  f  (t?VE  +  pt3VN)dt 

After  evaluating  the  right-hand  side  integrals  >f  Equations  (12),  (13),  (14), 
and  (15),  thi3  set  of  equations  can  be  written  in  the  matrix  form  as 
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TABLE  I.  VALUES  OF  ELEMENTS  OF  NORMAL  EQUATION  MATRICES 


Continuous  -  Time 

Least  -  Square 

Discrete  -  Time 

Least  -  Square 

%  Error  In 
Discrete  -  Time 

C2  =  4.608000xl06 

C2  = 

4.613271xl06 

.11 

C3  =  8. 294400xl06 

C3  = 

8.307429xl08 

.16 

C4  =  1.592525X1011 

C4  = 

1.595635X1011 

.20 

C5  =  3 . 185050xl013 

C5  = 

3 . 192629xl013 

.24 

C6  =  6. 552100xl015 

C6  = 

6.570238xl015 

.28 

atrix  in  Equation  (17)  be  denoted  G  and  its  elements  g„ .  Thus 

_  „  _  .3 

"  33  h 

3 

"  q21  =  g34  =  g43  =_^k 
4 

=  *  4  * 1,2  i  (i8) 

5  7 

=  g41  =  _g23  =  "g32  =  1J^k 

5 

=  q24  =  q31  =  ^42  =  0 

G'1  (19) 


gll  g22 


(20) 
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) 


Then,  the  elements  of  H,  denoted  h  ,  are  given  by 

hll  =  h3  =  ^22_ "  K1 
A 

h22  =  h44  =  gll  =  K2 
A 

i  .  ^12 

h12  =  h21  =  h34  =  h43  =  “  ~JT  =  ~K3  (21) 

h14  =  h41  =h23  =  'h32  =  ~q14  =  “K4 

A 

h13  =  h24  =  h31  =  h42  =  0 

where  K  ,  i«i  to  4,  is  defined  here  to  simplify  the  notation.  The  parame¬ 
ters  A1>1  i-i  to  4,  are  given  by 

(22) 

To  evaluate  Kf,  use  Equation  (18)  in  (20)  and  (21)  to  give 

K  =  2400  Q 

1  * - - 

❖ 

K2  =  2800 


*  ■  2', 


(23) 


After  steps  of  algebraic  simplification,  one  gets 


After  simplification, 


! 


(W2)k  =  (W2)k_!  •t'(vNftf  -  ‘  p(VEftf4  -  VEitl  ) 


4  k 


/VNf  "VNi^ 


12  'k 


+  u(\f  .  vEl\  (t4f  ♦  t?t£+  t2ft2.  *  tft3.  +  t4)k>  (29) 


20  /k 


Similarly  for  of  Equation  (16), 

(Vk  =  (Vk-1  +  (VEftf  '  VElt^  +  "(VNftf  '  VNlti) 

3  )k  V  4  /t 


~(VEf  -  VElj  <4  +  4*1  +  *f4  +  4\ 


12  k 


(VNf  ~  VNj|  ^f  +  tftl  +  tfti  +  tfti  +  tiJk. 


(30) 


Equations  (27)  to  (30)  are  recursive  updatinq  equations  for  Wj  which  are 
needed  in  Equation  (22)  for  the  recursive  determination  of  A  . 

Equations  (22),  (23),  (24),  (27),  (28),  (29),  and  (30)  constitute  the 
set  of  conti nuous-time  repression  algorithms. 


IV.  THE  SCRAPING  FEATURE 


Consider  a  regression  model 

V  =  Ak  +  Bk2,  k  =  0,  1,  2,  . . .  (31) 

where  V  is  measured  while  A  and  B  are  parameters  to  be  determined.  Assume 
that  B  is  much  smaller  than  A.  Under  this  condition  the  determination  of 
B  will  be  less  accurate  than  that  of  A  because  of  computation  round-off  er¬ 
rors.  ,0ne  way  to  get  around  this  difficulty  is  to  subtract  from  V  the  a- 
mount  Ak  where  A  is  the  a  priori  value  of  A.  A  is  related  to  A  by 

A  =  k  +  AA  (32) 

where  AA  is  the  error  in  k.  Thus,  a  new  regression  model  is  given  by 

V'  =  V  -  &k  =  AAk  +  Bk2.  (33) 

Usually  aA  is  small  enough  not  to  exert  a  dominating  effect  on  B.  For  a 
recursive  regression  algorithm,  the  part  of  Ak  known  a  priori  can  "scraped" 
off  from  V  recursively.  Therefore,  this  feature  will  be  called  the 
"scraping  feature".  A  demonstration  of  the  effectiveness  of  this  feature 
is  given  below. 


Example:  Consider  measurements  generated  by  the  model 
Vk  =  Ak  +  Bk2  +  nk,  k  =  1  to  100 

where  A  =  1.111111  and  B  =  .000001.  Measurement  noise  n^  has  a  normal  dis¬ 
tribution  of  mean  zero  and  variance  one.  Least  square  regression  with  and 
without  the  scraping  feature  are  used  to  estimate  A  and  B.  Results  are 
shown  in  Table  II.  Improved  results  due  to  the  scraping  feature  are  evident. 


Return  to  the  IMU  measurement  model  of  Equations  (1)  and  (2).  Here, 
accurate  determination  of  parameters  A£  and  A^  is  more  important  than  that 
of  Aj  and  A^.  The  scraping  feature  will  be  incorporated  into  the  regression 
algorithm  to  enhance  the  accuracy  of  the  estimated  and  A^.  Define 


VN  =  VN  "  M  =  AAll  +  A2t?  +  ,lA4t3 
VE  =  VE  "  =  +  A4t2  "  llA2t3 


(34) 
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TABLE  II.  LEAST  -  SQUARE  REGRESSIONS  WITH  and  WITHOUT  SCRAPING  fEATURE 


A  =  1.11111 
B  =  .OOOOOl 


Least  -  square  regression  with  scram'nq  feature 


K 

ESTIMATES 

f\ 

A  "  Error 

Error 

10 

1.1111  0.0000 

y.‘4r>8f)9xlf'-7  -5.4 

50 

1.11111  0.0000 

1.00060xl0'6  0.06 

100 

1.11111  0.0000 

0. 99795x 10-7  0.02 

Least  -  square  reqression  without  scraping  feature 


K 

ESTIMATES 

^  Error 

g  "  Error 

10 

1.11111  0.0000 

9. 53674x10” 7  -4.6 

50 

1.11112  0.0007 

8.34465xl0'7  -16.6 

100 

1.11113  0.0014 

8.34465xl0'7  -16.6 

where 


A 


1 


A 


3 


+  AAj 

^3  +  AA3 


(35) 


Equation  (34)  is  the  new  regression  model  for  IMU  self-alignment.  Compar¬ 
ing  Equation  (34)  to  (1)  and  (2),  one  notes  the  following  replacements: 


S 


A Aj  replaces  A^ 


AA^  replaces  A3 
Vj  replaces  VN 


replaces  V^. 

Therefore,  using  similar  replacements,  the  least-square  results  obtained  in 
Section  3  can  be  modified  for  the  model  of  Equation  (34),  From  Equation  (22), 
with  proper  replacement,  one  gets 


— 

- 

“  - 

3Ai 

K)  -K3  0  -K, 

a2 

as 

-K3  K2  k4  0 

W2 

AA3 

0  K4  Kj  -k3 

W3 

1 

< 

k 

-k4  0  -k3  k2 

k 

W4 

where  are  given  by  Equations  (23)  and  (24),  and  are  obtained  from 
Equation  '16)  as 

Wt  =  JtvN'dt 
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W2 

*  /(*;- 

t3VE')dt 

W3 

-  /tv-dt 

(37) 

W4 

■  * 

t3V')dt 

Once  (AAj)k  and  ( AA 3 ) k 

are  determined,  (Aj)k  and 

(A3)k  are  given  by 

iVk 

=  <Al>k- 

.1  +  (AAj)k  ’ 

> 

(38) 

<A3>k 

^A3^k- 

1  +  <M3>k  , 

The  algorithm  for  Wt  of  Equation  (37)  is  more  involved  than  that  given 
by  Equations  (27)  to  (30).  First  consider  the  scraping  of  VN  to  give  V^. 


(Vk  * 

tAl^k-ltk 

^Vk-1 

+  (AVk 

[<Vk- 

2  +  (AVk-J 

(tk-l  +  At) 

(Vk-l 

+  ^Al^k-2tk-l 

+  -  (A: 

- - ’ 

( ,AA1  )k-lfck 

( VN)k-l * 

(Vj-)^  can  be  obtained  from  (VE)k  in  a  similar  way.  Hence 

(VN}k  =  (VnVi  +  (AVk  ’  (Al}k-2At  -  (  Ml)k-ltk 

(VE}k  =  (V')^  +  (AVE)k  -  (A3)k_2/t  -  (AA3)k_1t|< 
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Next,  consider  the  effect  of  AA^,  and  AA3  on  W^.  From  Equation  (17),  re¬ 
placing  Aj  and  A3  by  AA^  and  AA3,  respectively. 


The  part  of  W  due  to  AA ^  and  AA3  alone  are,  therefore,  given  by 

<4Vk  ■  i  <aAi>k 

3 

(AW2)k  •  tj  (4Al)t  -  uj  (AA3)fc 

4  5 

<flU3>k  '  ik<4"3>k 

3 

(AH4)k  •  (AA3)k  *  ut®  (SA,)k 

4  5 

(AW  )k  are  the  a  posteriori  correction  for  (W^  after  (AA1)^,  (AAj)k,  and 
( AA 3 ) ^  are  estimated.  Ideally  AA1  =  AA3  =  0,  therefore,  the  corrected  (W^ 

denoted  (WC)k  a re  given  by 


(42) 
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(Vk 


(W3)k_i  -  (AW3)k_1 

'  T  (VEf  "  VeA  +  Vi  +  tl)k 
+  T  (VEf4f  '  VElVk*  (45) 

(w4)k_l  -  (AVk-l 

"  12  Af  "  VEi*k  A  +  tfti  +  Vi  +  ti^k 
+  7  AVf  -  ht'lh 

“iff  Af  "  Wk  +  tfti  +  tfti  +  Vi  +  ^1 


+  4  (VNf4  "  VNiti)k- 


The  complete  set  of  least-square  regression  algorithms  with  the  scrap¬ 
ing  feature  is  given  by  Equations  (36),  (38),  (39),  (41),  (43),  (44),  (45), 
and  (46). 

Testing  the  Scraping  Feature 


A  least-square  regression  alnorithm  with  scraoinn  feature  was  imple¬ 
mented  into  an  IMU  self-alignment  and  gyrocompass inq  software.  The  soft¬ 
ware  was  then  tested  on  a  long  ranne  missile  IMU.  A  total  of  1?23  sets  of 
incremental  velocities  AV^,  AV^,  and  A was  measured.  The  algorithm  was 
recursive  which  yielded  1222  sets  of  alignment  results  beninnina  with  the 
second  measurement  set.  For  the  sake  of  comparison,  another  test  run 
usine  a  least-square  alnorithm  without  the  scraoinn  feature  was  also  made. 

It  was  predicted  that  the  algorithm  with  scraping  would  result  in  fast¬ 
er  convergence  as  compared  to  the  algorithm  without  scraping.  This  pre- 


diction  was  confirmed  by  the  test  results,  although  the  improvement  is  small. 
Figure  1  shows  the  convergence  plot  of  values  for  «,  the  heading  angle,  D^, 
the  north  channel  drift,  and  Dp,  the  east  channel  drift. 

Algorithms  used  in  the  test  are  listed  in  the  appendices.  Appendix  A 
contains  the  algorithm  with  the  scraping  feature,  and  Appendix  B  contains 
the  one  without  the  feature. 

V.  CONCLUDING  REMARKS 

Two  software  features  have  been  described  in  this  report.  The  features 
are  applicable  to  self  alignment  and  gyrocompass ing.  The  first  feature 
is  the  use  of  contiguous-time  arithmetic  in  a  least-square  regression  algo¬ 
rithm.  This  featur*  .voids  the  error  caused  by  the  discrete-time  model. 

The  second  feature,  ca  ied  scraping,  is  intended  to  enhance  the  accuracy  of 
determining  the  drift  parameters  by  recursively  scraping  off  from  measured 
velocities  those  components  caused  by  platform  tilt  misalignment.  The  re¬ 
sult  of  scraping  H  to  reduce  the  effect  of  computer  round-off  error.  In 
the  sample  test,  the  improvement  of  gyrocompassing  by  scraping  is  small. 

This  is  due  to  the  small  round-off  error  to  begin  with.  In  the  case  where 
round-off  error  is  severe,  scraping  should  provide  appreciable  improvement 
in  gyrocompassing. 
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Figure  I.  Convergence  of 


APPENDIX  A.  FORTRAN  PROGRAM  FOR  LEAST  -  SQUARE  WITH  SCRAPING 


I  I  N4  ,t 


SUPklJUIlNL  AL13(lNt  ) 
t; 

C********** ****** ***************************** 

c 

C  s  lUKLD  At.  '‘4ALi3X‘'  UN  DISC 

C 

C**** *********************************** ****** 

c 

C . NLW  M1L0M  LEAST  SUUAKL  fit  ROUTINE  WITH  SCRAPITNO  t  LA  I  URL 

C  2  V  SLIT  80 

C 

C  UHL  kUUIINL  COMPUTES  LSI  THAT  Li-  Of  I1LI  AND  TILT  RAIL 

C  «Y  PLRI-OKM1NG  A  LEAST  SOUARES  I-  XT  OF  VELOCITY  DATA  Xfl 

C  KLf  LKENCL  LOOK  DINA  ILL 

C 
G 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCtCCUCCCCCCCCCCCCCCCCLCCCLLGL 

i: 

C  STANDARD  COMMON 

C 

C  7/30/74 

C 

C 

KLAL.  KiS.KU.Kl  ,  K  2  ,  K  3  ,  K  S  ,  K  I  .  KNOM  ,  KT  X  ,  LAMDA  ,  LAI  ,IDLNI  . 

*  MIRXA.MIR;<B,MIRXO,KTA,K ID.MUDIP 

C 

1NILGLR  GYPSL  .  A7GYP  ,  SLEW  ,  OD  H.C  ,  ZSt  Lti  .HLf  LG  ,CW  LG.EUt  G  ,Sf  LAV, , 

*  It:  M2  .  A  I-  LAG  ,  FS I  ,  YEAR  ,  DAY  ,  HOURS  .  SECS  ,  LRRCT  ,  AX  X  S 

*  ,  K  I LXO  ,  SYNC; ,  ERROR 
C 

COMMON  ANCLE  <  3 )  ,SMANG(6>  ,  K  1 S  ( 2 , 3 , 2 )  ,  Ttltl  <3>  ,CA(3)  ,  D  I  HtT  , GYPSL, 

P  A/GYP , SLLW<  6 ) , XMUUl , IMUU2 , £MUXN( 16) , XDXSC ,OD ISC ,DELC( 3)  , 

C  NUM(6>  ,  KDCL.  (  3  )  ,/t.f  LG  ,  HLf  LU  ,  CCt  LC  ,  MDf  LC.lkfLC.lt  f  LG  ,l.Rt  G<3)  , 

D  NECC  I  ,  I  Ht.A  <  3 )  ,K0(2 ,3) ,K1(2.3> ,K2(3> ,K3(3> ,K3(2, J)  ,KT<3> , 

L  KNUM<2> ,  K  I  X  (  3 ) ,DtLTA<6> ,DV(3) ,TURU(3> ,1DUM1 ,JUP ,SA(3>  ,L1DU. 

F  CCD  ,WL  (  3  )  ,ER<3>  ,OtLt  I  ,  R  DEL.  T  ,  CNT , RCN I ,  1RP(3 ,3)  ,  RMAX  , 

C  KMAXH, LAMDA, LAT (3) ,Dt  <3)  ,D1 (3)  ,D0(3) ,DG(3> (G,PtTA(6)  ,01 <3>  , 

H  J M X NR , IHLTI , XHET2 , DVR (3) ,  L.  I  D(  3 )  ,  ID0M2  ,  SFLAG  ,  OMEGE  ,  TRC  (  3 , 3 )  , 

J  TRD<3 , 3)  ,  I  PD  (3 ,3)  ,PH17  ,  PtllC  ,  DECK  A  ,  CONVt  .  1DUM3  ,AI  LAC  ,F  til  (  if)  , 

K  YtAR . MON  I H , DAY , HOURS , MINS, SECS , MSECS . COLD (2,3).  DCULD ( 2 , 3 )  . 

L  CNL.W  (2,3)  ,  DCNLW  (2,3)  ,0(2,3)  ,DC(2,3)  ,  LKKL  I  ,1ft  .MODIP.UIYIP  ,  It.  M2, 
M  JDLNI (3,3) , M I RXA ( 3 , 3 ) , M I RXB ( 3 . 3 ) , M TRXC ( 3 , 3) ,  IHD.IHF,  I  CP (3, 3) , 

N  I  CD  (3  .3)  ,  I  DR  ,K  IA(3)  ,KID(3)  ,  It  1UP(3,3>  ,  ItltTA(fc)  ,AX1S(3) , 

P  1  LSI  ,  1 1.  A  I  ,  FR  (  3  )  ,  I  HE  TM('j)  ,  1 1 D1  ( 2)  , 

UT 1 ML (200 ) . 1HNUR(2S0) , tHWES<2S0 ) ,E IDNU(2S0 ) ,t IDWS(2S0 ) , KK , 1 H3 ( 2S0 ) 
R  ,ALITP(12S) . KPLU I ,X1P(230> ,X2P(2S0> ,X3P(2S0> ,X4P(2S0) ,XSP(2SU> 
l,  ,  DNRCP  ( 1 20  )  ,  DNKC(J(  12S  )  ,  NSt.07  ,  ALP  I  0  (  12b  > 

I  .RTtlU  ,SYNC,NPU|- , IMUSN, ISPACE  , POF 1(30) ,BUF2(32)  ,  ERROR 
C 

cccccclclccc(.;cl(.lcccccccccccclcc(;lccccclccl(.cu.lccccclc(;ccccccccc(.(.cllc 

c 

DIMENSION  INLSt  (2) 

t  (JIJ1VALLNU  (  T  tit  I  M(  i  )  ,  XI  )  ,  (  I  til  I  M(  2  )  .  X2 )  ,  (  IHt  IM(3)  ,X3)  . 

A  (  I  HE  I M  (  4  )  ,X4)  ,  (  IHtlM(S)  ,XS)  ,  (  INLSt  ,  II.SI  )  ,  (INISF(2>  ,K  > 
MjIJIVALLNCt  ( t  R  (3  )  ,  t  KU  ) 
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DAI  A  L  1.5/it >77//H/ 

L  DtlLkMlNL  It  ALU  IS  IU  lit  LXLLUILD 

C 

it  (  INI  .Nt  .  U  )  UIJ  IU  1 
F  S  T 1 0  =  0 
KLIUK'N 
i: 

C  Dt  I LKH1NK  it  t  i.lU.1  UK  SLLUND  I  IML  IMk'UULU.  f  IKS  I 

C  I  HROULIt  WILL  NUI  LUM.'U  I L  I IIL  L.3F 

L 

1  Jt(tbli(llKM,4 
C 

C  F1K3I  I  IML  I  IlkUULtl  SLI  NUMBER  Ul  X  I LKA I  IONS 

L  AND  INI IXAL  CONDI f IONS 

L 

2  t SI  10  =  1 
NLSF  -=  1NLSF 

t 

C 

C  LIJNUER  I  fklJh  DLU/FIK  IU  RAD/SLL 

C 

t  F  RU=I  RU/( A600 . *DLCK A ) 

U=-FF  RU/3  . 

C 

I  1  =  0  . 

DTSQ=-DEL7  I  F  DLL  I  T 
D  I  LD=D  I  SU*I)LL  1  I 
DTUD=DTSQ*0f3G 
Ai  =  .  .?ii33S.4i33#D1SU 
A2=.25*DICB 
A,5=  .  2»D  !UD*U 
i: 

UNF=DUR< 1 ) 

VEF  =-DVR ( 2 ) 

VAF  =~DVR  <  vJ) 

C 

Ui=VNF  *Ai 
U2A=VNF*A2 
W2B-- VLF#AA 
UI3=UEF*A1 
W4A=UEF*A2 
W4£f=UNF*A.3 
US=VAF#A1 
K  =  l 

DX1  =  0  . 

0X3=0  . 

XI  =0  . 

XU  =  0  .  (I 
XJ  =0  . 

X33  =  0.0 
to  IU  s 

c 

C  SECOND  AND  SUBSLWULNI  UHLS  ttIKUULH,  (UMFU1L 

C  Ul  IU  US.  31  HI  33.  AND  XI  10  XS 

C 

4  K  =K  -t  1 


■SS3GS! 


I  1ME 
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r. 

I  IIIMI'UIL  Ui  III  OS 

c: 

1 1  1 1  mi  1 1 
i  f  - 1 1  miu 

l. 

ONI- ONI 

onf  =oni  mok-i  1 )  xii*i)iin  i»x i *  m 

OE I  -OIF 

Ot.l  OEI  I)0K<2)  XJiKbLlI  I  UX.4*  1 1 
OAI-OAF 

OAI -OAi  l)0k'<,4> 

C 

i isu- n  * i i 
TFCB=1|'3Q*IF 
1 1-  |JI>  -  II  SO*  1 1  so 
IF  01  1F30*(FCI< 

I  ISO-  11*11 
ricn=  I  130*1.1 
I  1  (J  0  -  1  n.U*  I  ISO 
c: 

I  2ND  1 1  SO -I  II  *1  It  I  ISO 
I 3K D-  I F  C B  + I FSO* I  1  *  I F*I 130  *  I ICO 

I  4  I  H-  IF  OUi  I  Fib*  IH  II  SO*  I  ISO)  II  *  I  IIH  I  I  UIi 

c: 

At-  .  3,43333333.4*  II  IK*DX1 
A 2  ■  I  h  011*0X1  - . 2*U*7  F  U I *DX  M 

A3  -  ,  3.433333.433*  II  (,H*DX3 
A4  =  .2S* IFQD*DX3+  .  2*U*1F0I *DXi 
C 

Ui=Wi  (  ONI  ON  1  )*  [2ND/6 .  I  (ONI  *11  SO  ON  1  *  I  ISO)*.*.,  At 
U2A  «W2A  -<  ONI-  -  ON  1  /  *  I  JR  0/ 1 «? .  M  ONF  *  I Fs H  ON  1  *  T  1 C  0  >/J. 

W2B=W2Eil  U*(  <  OIF  01.1  )*  I  4  I  li*  ,  (IS-  <  OF  I  *  1 1  00  01  1*1  100)*  .  2S  )  -  A,..' 
U.5=U.5  (Oil  OEI  1*  I2ND/6.  K  Oil-  *  11  SO -01 J  *  I  ISO)*  .  3-A3 
W4A=U4A  (OIF  Otl  )*.|  3KD/12 .  i  (Oil  *  1 1  C.lf-  01  1*1  U.IO/3  . 
W4B=U4b-IJ*<  (ONF  ONI  )*  I  4  IH*  .  US  -  (  ONF  *  1 1  0  0  ON  I.  *  I  1 01) )  *  .  IS  )  -  A4 
OS-  US  (OAI  0A1 )*13N0/6. I  (OAI  *IF SO  0A1*I  ISU)*.S 
C. 

C 

l  lUMF’U  I E  St  l(J  SS 

C 

l 

II  U-  1 F  #U 

1  FUSO  -  I  F  II*  I F  II 
PP  =  3S  ,  -•  64  .  *  I F  USO 

00=7. +s.*  rrusu 

KRKB  H  (.,1<*PP 
RkROR  R l< B*  IF 
RRR  l)=RRRl*  I  F 

C 

31=240  ,  *00/RRRH 
S2- 2H00  .  /RRR  D 
33=2100 . /RRKC 
S4=  16H0  .  * II /RRR B 
SS=  .5, /II  lH 

c: 
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U  LimiMIII:  Xi  10  Xb 

(, 

X iA-bi  *W1  B,3*W2A  64*W4A 

XI  B-  '3.5*U2b  04*W4B 

D  X  1  -  X 1 A  t  Xi  B 

X I I  XI 
XI  --XHDX1 

X2A=  -33#W1  4&2»W2A434#W3 
X2B--  32*W2B 

X2-=X2A+X2B 

X3A==  S4*U2A+Sl*U3-b,l*U4A 
X3B=  34*U2B  -33#U4B 

UX.i  -X.4A  t  X3B 
X33  ==  Xi 
X.i  ■ XJ+DX3 

X4A=-S4#W1  ~S3#U,i  BS2VU4A 
X4B~  02*W4B 

X4=X4A  t-X4£t 

xv=  s3*wb 

c; 

C  1  ILL  HUR  LK 

i; 

b  KILIU  =  12bB 
BUFK6)  -  13 
VU  10  l  -  l  .  3 
BUM (1+23)  -  LR(1> 

BUI  1<1  +  2<I>  =  MD(i> 

10  BUB  KH-6)  =  DVR  (  1 ) 

DU  20  1  --=  1,6 
20  BUF1  (  1  t  V  >  =  THt  I A  (  I.  ) 

DU  30  1  -  1,3 
30  BUB'l  < 1+iSi  =  I  HlilHd  ) 
BUM  (2/)  =  U  SB 
C 

JB  (K  . Lt.NLbl  )Rt.  I  URN 

r. 

ABLAU=1AND<  1 13, At  L A(, ) 

F3T10 =0 

IND 
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FORTRAN  PROGRAM  FOR  LEAST- SQUARE  WITHOUT  SCRAPING 
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■ULIl'U  I  I.NL  AL13dMI  ' 


L  4  *  4  4  4  4  4  *  4  X. 4  * 4  >i.  X  4  4  *  4  *  *  *  V  4  4  (  4  *  4  *  4  *  4  4  >1  4  4  V  4  *  *  &  X- 4 4 

c 

C  .  1UKI  D  Ail  "CALiXf."  UN  U  UX 

L  KLV.  20  0 L  P  !  Lit) 

L 

L44444444444444444444444.4C  X  4444  44444C  444444444 
I. 

C  NOW  M1LUM  LLAOI  OUUaKL  ill  KOOllNt 

I.  /AUGAi  ROY.  1.JAUGU0  i. 

C 

l  UllO  KUUl  INL  CUMI'OIIX  tXliMAlli,  Ul  I  j  1. 1  AMI)  111.1  MAIL 

C  BY  PLNS-ORMINU  A  LLiV.il  SMOARlX  i  tl  III  L'.X  CJLI  L  Y  Oh  I A  IN 

(,  KLILKLNLX  LOOK DINA  1 1  1» 

L 

C 

C.XXL(XXLCXU.XLXCXLLXLLLLLLL(XLUXL(XC.lXLLLLl.LLLLLLCXLL(XL  LlXtXUX'.XUX 
C 

C  0  I  AN  DAK  I)  COMMON 

L 

L  s/30/74 

t: 

L 

KIAL  K  1 0  ,  K  0  ,  K  1  .  K  2  ,  K  .  • .  K  0  ,  K  I  ,  K  NUli ,  K  I  X  ,  I  AMD  A  .  LA  I  ,  1DLNI  . 
t  M  1 R  XA  .  H  I  KXH  .  rt  I  K  Xi .  ,  K  I  A  .  K  1 1)  .  HOD  1 1’ 

C 

i  M  I  LC>l.  K  OYI'GL  ,A/OYI'  M  LW  .  01)100  . /XI  LI.,, Ill  I  I  l,, CXI  IX. LKI  U.CI  L.  AC, , 

%  ICM2  ,  Al  L  MU  ,10  I  .  I'LAK  .  DAY  ,  IIOURO  OtXS.LHRI  I.AXIX 
l  ,K’U  10  XYNL  .LKKOK 
C 

COMMON  A  Nl>  I  L(3>  XMAN0(6>  ,K10<2  .>.2)  I  UL  I  <3>  XA(3>  .  D  I  ML  i  XYLOL  , 

B  A/XYL  ,  OL.LU  (  6  )  ,  mom  ,  IMUU2  .  JrtOlNC  16)  .  11)1  OC  ODlOL  .  I)tLC(  ,J  )  , 

(.  NUMC  A)  ,  KDLL  (  3)  ,  /XI  I  L  ,01  I  10  XU  L6.M0I  LX  ,  1KI  LX.  II  I  LX  ,LKI  0(3  )  , 

D  NEGU  I  .  I  HI.  AC  3)  ,  KU  (  .  .5 )  .  K  1  (  J  .  3  i  ,  K3(3  >  ,KOL'  .3)  .  K  I  (  3  >  , 

L  KN0M(2>  ,  K  I  X  (  3  )  .DLL  I  A  (  6  )  ,  DO  <  3  )  ,  I  UK  0(3)  ,  lDIIMl  ,  JUC  ,‘.A(3>  ,  L.  I  DO, 

f  CCD  , UL ( 3 )  ,;_R<3>  .01  I  I  I  KDLL  I  .  t.N  I  ,KCNI  ,  I  l(l>  (  3  .  ..<  i  ,  K  MAX  . 

C  K  MAXI.  .L.AMDA  ,LAI  ( 3 )  Dl  C3>  ,DJ  (3)  ,1)0(3)  ,D0(3)  .  C,  .HI  IA(6>  XI  (3)  , 

0  I  MINK  ,  I  III  1 1  .  I  ML  .12  .  DVK  <>)  XI  !)<•$)  .  1 DUM2  ,  Cl-  LAC  ,  UMLGC  ,  IRC  <3, 3  >  . 

J  IK  l)(  3 , 3)  ,  7LD<  3 .3)  ,1'IIJ  i.  .I'll  U  ,1)1  OKA  ,  LUNCH  .  iDLIM3.nl  LAO  ,101(13), 
l<  Yl.AK  ,  MON  III  ,  Day  .11001(0  ,  Ml.  NO  ,  OLLO  .  MOLLO  .001  0(0.3!  ,  DC.ULIHO  ,  J  )  . 

(  LNLW  ( DLNl  W  (  2  ,  .■  >  ,0(2,3)  .00(2,3)  .  LKKL  I  ,  I.  t'l  ,  HOD  It’  ,  U  I  Y  II'  .  1 1  M2  , 

M  X  DON  !  (  3 , 3  )  ,  M  I  8  XA  (  S  ,  5  )  ,  Hi  «  XB  (  3 , 3  )  .  M  1 1<  aC.  (  3  ,  S  )  ,  I  HD  ,  1  HI  .1  LB  (33), 

N  I  0D(  ,J  ,  3)  ,  I  DK  ,K  1  A(3  <  ,K  I  Di  3  )  .  1 1  101'  (3 , 3)  ,  I  UL  1  A(0  )  ,  AX10(  3)  , 

!>  ILOI  ,  I  LA  I  ,  I  R  (  5  )  ,  IIIL.  IMCj)  ,1.  I  DM2)  , 

01  IHL(L'MI)  ,  IHN0KC20U)  .  I  HULL,  (  20 II  )  ,t  I  UNO  (  200  )  ,LU>WO(20l)>  ,  KK  ,111.,  (20(1) 
H  .  ALL  f  H  (  1,30)  ,  KI’L.I)  I  ,  XIIM20U  )  ,  X2I>  ( 200 )  .  X,JL  (20  0  )  ,X4P(2\»U  )  ,  XSL(2'UI  ) 

0  ,  DNKC.I'  ( 1 20  )  ,DNKLW(  120)  ,  NOLI)'/  ,  ALL  10(1 20 ) 

r  ,R  U.IIJ  ,'OYNC  ,NBUI-  .  1M00N  .  101'ALl  ,  BOF  1(30)  .Bill  2(3,3)  .fRROH 

(, 

LXLXLLLLLLLLLL l.l.X(  LLLLLI  l.XLUXXl  C,f  I.CXL.LU  L.LXLIXLLXLU. (,(..(  XLXUXLXLLLUXL, 

D 1  ML  NO  KIN  J  ML 01  (2) 

lOUlVAILNIl-  (  Hit  IM<  )  )  .Xi  )  (  I  l!L  I IH  2  )  .  X2  )  ,  (  MIL  I  IK  3 )  X3) 

A  (  IIIL  I  M(  4;  .  X4)  ,  (  I  HI  Il1<  o  i  .  X.)  (  l  Nl.iil  ,  I  LOi  )  ,  I  INL.Oi  i  2)  ,K  ) 
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LUlUVNI.CNC;  <i  IMS)  .i  KM) 

daia  ii4/u-/’/y/b/ 

t.  I)L  I  E.KM1NL  II  Al  1 S  I'j  Ml  Ht  LX  •MU  1 1  D 

(. 

J I  (  INI  .  ML  .  U  )  Ml  I  U  J 

i  s  I  i  u  n 
Kt  nii.N 

L 

C  DEI  LKM1NL  it  t  JkLI  UK  U  LUND  1  i  ht  HIKUUUi;.  tlKOI 

C  litRUOUU  UII  L  NUI  COMi'IITE  1111  I.SE 

C 

)  ltd  Mil)  j  >  2 , 4 , 4 

c: 

L  t  i  K 1. 1  11ML  I  tlKUUUtl  !>l  I  NUMUI  k  IJI  iU.kAUONO 

C  AND  INI  IJAI  CONDI  I  i  ONL> 

L 

2  I  Lllll  -i 

NLCil-  -  l  Ni  !:if 
L 
C 

t:  UINVL.k  I  t  KUM  DLU/fik  I U  RAD/I.L  C 

C 

I  FKU-  I  KIJ/<S60  0  .  #Dt(,KA> 

r: 

C  M  l  MlJt#  LION  INOLRUD.  MAY  Id.  Kt.  MOULD  Will!  iON .  VX 

c 

U  -  ttkU/S, 

i: 

l  i  =  o . 

DTSU=DEL  I  I  *DELH 
DTCli-DT  OUXDLL  I  t 

d  ruo- 1)  rsci*i>  tub 

Al  I  iM 

,  ;_".j*D7 CD 
AS-  .  L'#D  I  UD*U 
C. 

VNF  =  DVk(  1  ) 

VEFr*  DVR  <  2 ) 

VAt  ~  DVK  ( S  ) 

C 

W1  ==VNI  *A1 
W2A«VNF  #A2 
WLba-  VLF  *AS 
UJMVEF*A1 
W4A- Vtt #A2 
U4H-VNI  *A3 
WS--VAI  *A1 
K  =  i 

UU  TO  «, 

l. 

SECOND  AND  MJDCF  (JIILNI  I  I  Ml  0  IHKOUUH,  COMIUlt 
C  W1  10  US  01  10  S»b.  AND  XI  III  X5 

C 

' !  K--Ki  1 
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l  J  ML. 


L 

L 


CUMHIlt  Wi  10  WS 


I  L  :  1)1  l  I  I 

I I  i DLL  I  ! 


VNL-- VNI 

VNI  =  VN)  H)V.((  1  ) 

'.'l  l==Vl  F 

VES'-Vti  DVR (2) 

VAl  =  VAl 

VAI  -VA1  DVR  <  J  > 

I  F  30  1 1  *  II 
i  h c &  =-  i':-ho*if 
I  (  01)  If  SUH.  1 1  30 
I  IS0  =  11*11. 

I  i  ( ,  It  I  130*  I  I 
I  101)=-  I  F30*  I  ISO 

l,.'ND  II  30:  I)  *1  III  1  SO 

r.3l<D  (I  CH-t  IF  30*1  F+1F » 1130; F lLfi 

I  4  I  H-  IF  01)1  IF  ( ; l: *  I  1  i  I F  30*  I  130;  I F  *  I  Il.iH  I  I  01) 

Wl-Wl-  (VNF  •  VNI  )  *  I  L'ND/6  .  ■*  (  VNF  *  IF  30  VN)  *  I  130>* . 3 
UiiA- U.-’A  (VNI  VNI  )*  F.5RD/12.  +  (VNI-*TF  LI)  -VNI *1  11.11 )  /3  . 
WL'H-Ucl: : IJ*  (  (  VLF  VI.  I  >*  I  4111*  .  03  (VI.  F  *  IF  01)  VL  I*  I  101))*  .21 
WJ--W.S  (VI. I  -VL.  I  )*lc'ND/6.  K  VLF  *11  30  VE  1*11  30 )  *  .  3 
U4A  U4A-  (VFF  VI  1  )*  I  JKD/IL! . -;  (  VU  *  IF  l.H  VI.  I  *  I  1L1)  )/,'>. 
W4B-U4K  II*  (  <  VNI  VN  l )  *  F  4  T 1 1*  .  03  •  (  VNI’*  1 1  01)  VN I  *  F  IUl)  >  *  .  1?' 
W3-W3-  (VAI  VAX  )*i;-Nl)/(,  .-!  (VAF  *  I F  30  VAl*II30)*.3 


LllhF’Ull  3 ;i  III  33 


l F  U- 1 1  *U 
1FU30  II  LI*  It  U 
KF’-33  .  -I  64  .  *  I F  1)30 
00  =  '/.  13.*  OHIO 
kkkt:  I F  LI:*I'F' 
kRRORRRK*  IF 
kkki)  kkkL*  if 

31=24(1  .  *00/KKkL: 

S2=2l)0ll.  /RkRD 
33-21  00  .  /  kkkL. 

34  =  1690  .  *IJ/kkkfi 
33-  3./ II  LH 

(OWU  It  X 1  111  X3 

X 1  A  =  3 1  *W  1  33*u;.  A  34*W4A 
X1H=  3.1*021)  3  4  ♦  W  4  F< 

X1  =  X1AIX1I: 

X2A---3  J*WJ  *32*0, JA  i  ',.4*0,1 
XL'k-  3L'*W;.'li 

x2=-x.ia-*x;?& 


X3A=  U4<U2A*ai*W.5 -S.5*W4A 
X3B~  U4*U2B  ‘,3*W4B 

X3=X3  A+X3B 

X4A--  U4*U1  33*U3  ( i>2*W4A 
X4B-  U2*W4B 

X4-X4A  •  X4B 
X'5*  3'»*U'd 
L 

l.  (ILL  BUI  F  l  K 

c 

XltiU  =  123D 
BUF 1(6)  =  13 
DU  i  0  1  1,3 

BUF  i  (  I  +23  >  =  ER(I> 

BUI  1(1+20)  =  l  ll>(  1  ) 

1(1  BUF  1(1+6)  -  DVR  (  1  ) 

DU  20  1  •-  1.6 
20  BUF  1  <  1+9)  =  7HI-  I  A(  X  ) 

DU  30  1  1.3 

3(1  BUF  1(1  +  15)  =  TMi  TM<  I  > 

BUF 1 ( 2/ )  -  ILUF 

i: 

IF  <  K . LL.NLUI  >KtIUKN 
C 

AF  LAU=1AND(113,AFLAU) 

F  STi  0-=0 
LND 
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